Impact of a Warm Dark Matter late-time velocity dispersion on large-scale structures 
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We investigate whether the late-time (at z < 100) velocity dispersion expected in Warm Dark 
Matter scenarios could have some effect on the cosmic web (i.e., outside of virialized halos). We 
consider effective hydrodynamical equations, with a pressurelike term that agrees at the linear level 
with the analysis of the Vlasov equation. Then, using analytical methods, based on perturbative 
expansions and the spherical dynamics, we investigate the impact of this term for a IkeV dark matter 
particle. We find that the late-time velocity dispersion has a negligible effect on the power spectrum 
on perturbative scales and on the halo mass function. However, it has a significant impact on the 
probability distribution function of the density contrast at z ~ 3 on scales smaller than 0.1fa _1 Mpc, 
which correspond to Lyman-Q clouds. Finally, we note that numerical simulations should start at 
Zi > 100 rather than %i < 50 to avoid underestimating gravitational clustering at low redshifts. 



PACS numbers: 98.80.-k,95.35.+d 
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I. INTRODUCTION 

In the standard cosmological A cold dark matter 
(ACDM) scenario, most of the matter content of the Uni- 
verse is made of CDM particles, which are cold and col- 
lisionless. This means that they have a negligible veloc- 
ity dispersion during the matter-dominated era and den- 
sity fluctuations on almost all scales (except very small 
scales at early times) grow through gravitational insta- 
bility. This leads to a hierarchical scenario for the forma- 
tion of large-scale structures, as the amplitude of density 
fluctuations at early time (e.g., at the beginning of the 
matter-dominated era) is larger on smaller scales. Then, 
small scales turn nonlinear first and merge to build in- 
creasingly large and massive structures as larger scales 
become nonlinear in the course of time [l[ . This scenario 
(with an extra dark energy component or cosmological 
constant) is in good agreement with a large variety of 
cosmological observations, such as the cosmic microwave 
background (CMB) and galaxy surveys 

However, this CDM model may disagree from observa- 
tions on small scales (below the size of galaxies). Thus, 
CDM simulations typically predict too many satellite 
galaxies around Milky- Way-sized central galaxies as com- 
pared with observations [J-Q- They also predict power- 
law density profiles, p ~ in the center of virialized 
halos p} , whereas dark-matter-dominated dwarf galaxies 
[H and some disk galaxies @ exhibit flat density cores. 
This is the so-called "core-cusp problem" [Tc| ■ 

One possible solution to these small-scale problems is 
a warm dark matter (WDM) scenario, with dark matter 
particles of a mass on the order of IkeV. This intermedi- 
ate case between the "cold" and "hot" dark matter sce- 
narios provides a non-negligible velocity dispersion and 
a significant free-streaming that erases density fluctua- 
tions on small scales (mostly during the period where 
the particles are relativistic) . This helps to cure the 
small-scale problems of the CDM scenario, while being 
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This fa- 
vors a mass on the order of IkeV [15j, [lfj. For larger 
masses we recover the CDM scenario and for smaller 
masses we recover the hot dark matter scenario, where 
structure formation begins too late (in particular, this 
is ruled out by the Gunn-Peterson bound [13]: quasar 
spectra show that the Universe must have been reionized 
before z ~ 6, which requires galaxy formation by this 
time). Similar lower bounds on m are also obtained from 
the observed velocity dispersion of dwarfs galaxies and 
from the Lyman-cv forest [l8l - [2l| . 

We must note that these small-scale problems may 
also be cured by the physics of the baryonic compo- 
nent, within the CDM scenario. For instance, reion- 
ization of the intergalactic medium [2^, [H| or feedback 
from stars and supernovae (2~4| suppress star formation 
in small satellite halos. Only a small fraction of the low- 
mass dark matter satellites would then shine in the sky 
and appear in galaxy surveys. This would reconcile the 
observed abundance with the CDM prediction but there 
remains some discrepancy for the shape of the satellite 
luminosity function [25]. Supernovae explosions may also 
transform a cusp density profile into a cored one, within 
small dark matter halos [26H28I ]. However, it is a difficult 
task to check that such models can explain galaxy prop- 
erties from massive to dwarf galaxies and from z = to 
higher redshifts [H. 

Therefore, WDM scenarios remain interesting alterna- 
tives to CDM that are still being investigated in many 
works. As recalled above, a particle mass on the order 
of 1 — lOkeV is a good candidate and it may correspond, 
for instance, to sterile neutrinos [30l - l36| or to gravitinos 
[37], [38| ■ At early times and on large scales, the formation 
of large-scale structures within WDM scenarios is stud- 
ied through the linearized Vlasov equation [39l - l43j . At 
low redshift and on small scales, the nonlinear regime of 
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gravitational clustering is investigated through numerical 
simulations frll HI] and halo models based on such sim- 
ulations j4^,|4g|. In practice, one often uses the same N- 
body codes as for CDM scenarios and the only difference 
comes from the density power spectrum that is set at the 
initial redshift Zi of the simulations. This means that one 
takes into account the high-fc cutoff due to free-streaming 
during the relativistic era but neglects the nonzero veloc- 
ity dispersion at low redshifts, z < Z{. This is legitimate 
because the relative importance of this late-time velocity 
dispersion decreases with time and the main difference 
between the CDM and WDM scenarios with respect to 
large-scale structures arises from the high-Zc cutoff of the 
power spectrum built during the relativistic regime. Nev- 
ertheless, it would be interesting to have a quantitative 
check of this approximation. This is the goal of this pa- 
per, where we obtain a quantitative estimate of the im- 
pact of the late-time WDM velocity dispersion on the 
formation of large-scale structures. We also estimate the 
sensitivity of the gravitational clustering measured at low 
redshift on the initial redshift Zi of the simulations. 

Here we do not consider the inner regions of virial- 
ized halos, where the finite velocity dispersion can have 
important effects because of Liouville theorem. Indeed, 
this implies an upper bound on the coarse-grained phase- 
space distribution function [13], which can lead to cored 
density profiles instead of cusps 48] (but the behavior 
in central regions remains difficult to predict (49[). In 
contrast, we consider the cosmic web, that is, moderate 
density fluctuations or large scales, as well as the halo 
mass function itself. Then, using perturbative methods 
or the spherical dynamics, we compare such statistics 
(the power spectrum on large perturbative scales, the 
halo mass function, and the probability distribution of 
the density contrast) between the CDM and WDM sce- 
narios, where we neglect or take into account the late- 
time WDM velocity dispersion. To simplify the analy- 
sis and to go beyond the linear regime, we use effective 
equations of motion similar to standard hydrodynamics. 
They involve a simplified pressurelike term in the Eulcr 
equation, associated with the late-time velocity disper- 
sion, that is chosen so as to agree with results from the 
Vlasov equation at linear order. This should be suffi- 
cient for our goal, which is only to estimate the order of 
magnitude of the impact of this late-time WDM velocity 
dispersion. Thus, our study is complementary to Ref.plj 
who investigates the effects of the velocity dispersion at 
low redshift through the linearized Vlasov equation. Our 
approach is not exact, because we use a fluid approx- 
imation, but it allows us to consider nonlinear density 
fluctuations. In particular, our goal is not to study in 
accurate details a specific WDM model but to investi- 
gate the generic impact of a late-time velocity dispersion 
on the cosmic web. After describing these effective equa- 
tions of motion and our approach in Sec. [TTl we present 
our results for a IkeV dark matter particle in Sec. Mil We 
also consider the impact of the choice of initial redshift 
in numerical simulations in Sec. IIVI and we conclude in 
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II. EQUATIONS OF MOTION 

A. Effective Euler equation 

During the matter-dominated era, on scales much 
smaller than the Hubble radius and for nonrelativistic 
particles, the evolution of the dark matter density per- 
turbations is governed by the nonrelativistic Boltzmann- 
Vlasov equation [l[, 
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Here r is the conformal time, a(t) the scale-factor, x and 
p are comoving coordinate and momentum, m is the mass 
of the particles, n the mean dark matter number density, 
and /(x, p, t) the phase-space distribution function. Be- 
cause the gravitational potential $ depends on the dark 
matter density fluctuations (here we neglected the con- 
tribution from baryons, or, more precisely, we neglect the 
difference between baryons and dark matter), the Vlasov 
equation (|TJ) is nonlinear over /. Then, one usually lin- 
earizes the system ([I])-© to study the evolution of the 
dark matter perturbations [3^, H3] ■ 

The phase-space distribution /(x, p, r) includes the 
distribution of velocities of the dark matter particles at 
any point in space. However, for observational purposes 
one is mostly interested in the density and peculiar ve- 
locity fields, which are the lowest-order moments of /, 
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By taking successive moments over p of the Vlasov equa- 
tion (p} one obtains an infinite hierarchy of equations 
that no longer depend on the velocity coordinate [l|, [5(| ■ 
However, contrary to the case of collisional fluids, this 
hierarchy cannot be closed by writing the second-order 
moment (the velocity dispersion) as a pressure that obeys 
a definite equation of state. In the case of cold dark mat- 
ter, this hierarchy can still be closed by neglecting the 
velocity dispersion (whence the name "cold" ) and one is 
left with only two fields, the density p(x, r) and the mean 
peculiar velocity v(x, r). In the case of warm dark mat- 
ter, the velocity dispersion cannot be neglected at high 
redshift and one must work with the Vlasov equation. 
This large velocity dispersion leads to a significant free 
streaming that damps the early time density power spec- 
trum at high fc, as compared with CDM, because parti- 
cles can escape from small potential wells [ll|. This is 
very efficient at early times, when the particles are still 
relativistic. Therefore, the damping is greater for par- 
ticles of smaller mass m, which become nonrelativistic 
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later. Afterwards, because the typical comoving veloc- 
ity dispersion decreases with time because of the Hubble 
expansion, these effects become less important on large 
scales. (However, because of Liouville theorem, the fi- 
nite velocity dispersion, which leads to a finite value for 
the initial phase-space density / -as opposed to the cold 
case where it would be a Dirac distribution- leads to an 
upper bound for the coarse-grained phase-space distribu- 
tion [47j • This can have a significant effect on the matter 
distribution on small scales, within collapsed halos, even 
today HI.) 

To go beyond the linear regime, one usually runs nu- 
merical simulations to investigate the formation of the 
cosmic web and collapsed halos. In practice, one often 
uses the same codes as for CDM [n|, Q M, EH and 
the only difference is encapsulated in the initial density 
power spectrum set at the initial redshift Zi of the simula- 
tion (i.e., the high-fc cutoff due to early time free stream- 
ing). Sometimes, in order to mimic the effect of the up- 
per bound for the coarse-grained phase-space distribu- 
tion, one adds an initial white-noise velocity component 
to the particles of the simulation [52h54| . However, this 
is not truly equivalent to the full phase-space distribution 
/(x, p, r) because the "particles" used in the simulations 
are actually "macroparticles" , which model large clumps 
of matter, rather than the physical dark matter particles. 

Our goal in this paper is to estimate the impact of a 
late-time WDM velocity dispersion on large-scale grav- 
itational clustering. As explained above, within a hy- 
drodynamical approach, one cannot truly close the hi- 
erarchy of equations obeyed by the successive moments 
of the phase-space distribution. However, because the 
fluid equations are much easier to follow beyond the lin- 
ear level than the Vlasov equation and we only look for 
an order-of-magnitude estimate in this simple study, we 
will use a simple closure inspired from linear theory. As 
shown for instance in [39l |40| , in the matter dominated 
era, one can derive from the linearized Vlasov equation 
an evolution equation for the matter density contrast of 
the form 



d 2 6 
dr 2 
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-n m H 2 - k 2 c 2 



S[S,t] 



(4) 



where 51 m (r) is the matter cosmological parameter, % — 
a is the conformal expansion rate, and we denote with 
a tilde Fourier-space fields such as #(k, r). This corre- 
sponds to Eq.(2.85) in (39[, with a different time variable. 
Equation Q becomes identical to the standard equation 
for the linear density modes in CDM scenarios when we 
set c s = and 5 = 0. 

The new term in the left-hand side is similar to a pres- 
sure term (at the linear level), and c s would be the co- 
moving sound speed. The comoving Jeans wave number 
kj would thus correspond to the comoving free-streaming 
wave number fcf s , with 
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However, this is only a formal analogy because there 
is no true pressure as we consider a collisionless fluid. 
This term arises from the nonzero velocity dispersion 
and its evolution with time is set by the comoving free- 
streaming wave number (J5]) rather than by a thermody- 
namical equation of state. At this linear level, kf s (r) is 
obtained from the analysis of the linearized Vlasov equa- 
tion. This yields [H, 5<| 



fcf s (r) = fcf s (r cq )./ , 



(6) 



where et oq is the scale factor at the matter-radiation 



equality, a cq = 1/(1 + z cq ), with z c 



kfsijeq) 
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Mpc" 1 , 



(7) 



where y 2 ~ 10 depends on the shape of the initial veloc- 
ity distribution, g& is the effective number of relativistic 
degrees of freedom at decoupling, and m is the mass of 
the dark matter particles. 

The new source term in the right-hand side of Eq.Q 
can be decomposed as a memory term Snb (i-e., an in- 
tegral over past times) and an inhomogeneous term Sb 
(i.e., an integral over the "initial" condition at matter- 
radiation equality) ; see [lij H(| ■ The nonlocal term Snb 
decreases as (k/kf s ) 4 at low k, hence, it is subdominant 
on large scales as compared with the new term in the 
left-hand side. The inhomogeneous term only decreases 
as k 2 , but its relative importance decays with time, as 
compared with the growing mode associated with the 
left-hand side. 

In this paper, we neglect the source term S, and we 
investigate the equations of motion 

|*+V.[(l+*)v] = 0, (8) 

OT 

f +Hv + (vV)v = -V$-^, (9) 

OT p 

where the gravitational potential is given by the Poisson 
equation ([2]), which also reads as 



v 2 $ = -n m n 2 5. 



(10) 



(5) 



Equation (jSJ is the usual continuity equation, which is 
exact when v is the mean peculiar velocity as defined by 
Eq.([3|). Equation (|9|) is the Euler equation, which only 
differs from the CDM case by the last term. This equa- 
tion is not exact, because as explained above this last 
term should be written in terms of the nonzero velocity 
dispersion, which obeys a third evolution equation that 
involves the third-order velocity moment, and so on. The 
expression ([9]) is the simplest closure of this hierarchy 
that agrees with the linear theory ^ (where we neglect 
the source term S). In a fluid analogy, where c s would 
be a uniform isothermal sound speed, this term would 
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be modified beyond linear order because it would read as 
{Wp)/p instead of (Vp)/p. However, because this anal- 
ogy is not exact (and does not need to be strictly fol- 
lowed) and to simplify the analysis, we keep the linear 
term of Eq.Q. This also ensures that we always mimic 
the slowdown of gravitational collapse on small scales as 
compared with CDM. (In contrast, a truncation at a fi- 
nite order over 5 of V<5/ (1 + 5) is not very well behaved. 
For instance, depending on whether the truncation order 
p is odd or even a term (— l) p ~ 1 8 p ~ 1 V(S either slows down 
or speeds up the dynamics of overdense regions). 



B. Setup 

The closure ([9]) should be sufficient for our purpose, 
which is not to obtain the most accurate predictions 
for large-scale statistics but to estimate the impact of 
a late-time velocity dispersion. Equations ©-(O read in 
Fourier space as 



with i/ — 1.12 and 



as 

8t 



(k,T) + 0(k,r) = -/ dk!dk 2 <5 D (k a +k 2 -k) 

xa(k 1 ,k 2 )0(k 1 ,r)5(k 2 ,r), (11) 



°— (k, r) + W(k, t) + [l + e(k, r)]6(k, r) 



dkidk 2 fe(k x + k 2 -kJ^ki.kaJ^ki.r^ka.T), 

(12) 



where we introduced the velocity divergence, 9 — V • v. 
The usual kernels a and /3 are given by 

(k 1+ k 2 ).k! Ikx+kafQci-ka) 
a(ki,k 2 ) = ,/?(ki,k 2 ) = - 



and the factor e(k, r) writes as 



e(k,r)=- 



2A^2 ^2 



fcfs(r) 2 



(13) 



(14) 



Combining Eqs.lfTTj) and fj 12[) at linear order, we recover 
the left-hand side of Eq.Q. The equations of motion 
(fTTj) - (ll2p have the same form as those studied in 55 1 
in the context of modified gravity models, but with a 
different kernel e(fc,r). Therefore, we can use the same 
methods as in [55j to investigate large-scale structure for- 
mation. 

To fully define our system, we must specify our initial 
conditions. Following [Hjj . we write the linear density 
power spectrum today in terms of the reference CDM 
power spectrum as 

P L ,WDM{k 7 z = 0)= P L ,cou(k, z = 0) T(fc) 2 , (15) 

with 

T(k) = [l + (afc) 2 T 5/ ", (16) 



a = 0.049 



f— V 

VlkeV/ 



0,. 



0.25 



h 

0T7 



Mpc. 
(17) 

The high-fc cutoff (fTB]) is due to free streaming in the rel- 
ativistic regime, which damps the growth of small-scale 
clustering as compared with CDM. 

Then, in a manner similar to what is done in numeri- 
cal simulations, we choose an "initial" redshift Z{ = 100, 
and we write the initial power spectrum as Pt{k,Zi} = 
[D + { Zi )/ D + {Q)]P WDM (k,0), where D + (z) is the CDM 
linear growing mode, which is independent of k (it is the 
growing solution of Eq.Q with c s = and S = 0). Next, 
for a given initial redshift Zi, we compute the subsequent 
formation of large-scale structure defined by the equa- 
tions of motion (ITl"1) - (fl"2")) with cither e ^ (i.e., taking 
into account the late-time velocity dispersion through the 
last term in Eq.©) or e = (i.e., neglecting this late- 
time velocity dispersion), starting with the same initial 
WDM power spectrum (fl"5|) . Comparing these two re- 
sults we obtain an estimate of the impact of the nonzero 
velocity dispersion at low redshifts, below Zj. Comparing 
with a pure CDM scenario, with the initial power spec- 
trum Pcdm, we compare this effect with the damping due 
to the initial high-fc cutoff (|16p . which mostly arises from 
the relativistic era. This will allow us to check whether 
the late-time dark matter velocity dispersion can be ne- 
glected on large scales (outside of collapsed objects where 
the upper bound associated with Liouville theorem may 
be important). This is often the assumption used in nu- 
merical simulations, where the system is modeled with 
the same codes as for CDM and the only difference lies 
in the initial condition at redshift Zi , that is, in the damp- 
ing (Unj) of the initial power spectrum [IJ |4J, |46|, |51( . 

In some cases, one adds to the macroparticles used in 
the N-body simulations an additional initial random ve- 
locity, f rms , drawn from the thermal distribution f(v) of 
the WDM particles [52M541 ] . However, this is not fully le- 
gitimate because these macroparticles have a much larger 
mass (~ 10 5 M Q ) than the WDM particles and a clump 
of this mass of many elementary dark matter particles 
would have a much smaller (almost zero) mean velocity. 
In fact, setting such initial conditions at Zi = 100 leads 
to spurious power at high k in the power spectrum mea- 
sured at later times (starting at Zi = 40 _appears to avoid 
this problem because u rms is smaller) [521 ] . This can be 
understood from the fact that adding these random ve- 
locities is equivalent to model a CDM scenario, with the 
damped power spectrum (|15|) . to which is added a k 2 
high-fc tail associated with the random velocity compo- 
nent. Indeed, this adds a white-noise component to the 
initial velocity power spectrum, which corresponds to a 
k 2 tail for the density power spectrum (when we decom- 
pose over growing and decaying modes). On the other 
hand, these additional random velocities can be seen as 
an effective tool to build an upper bound on the phase- 
space distribution function and to investigate the forma- 
tion of cored profiles associated with Liouville theorem 



5 



|5q . However, this shows that it is not easy to include 
the WDM velocity dispersion in numerical simulations 
in a realistic fashion. This is another motivation for the 
analytic study presented in this paper. 

In the following, as in [44j , for numerical computations 
we adopt a background cosmology that is consistent with 
WMAP7 H, Q m = 0.2726, Ct A = 0.7274, il h = 0.046, 
h = 0.704, n s = 0.963, and a s = 0.809. 

We focus on the case of a 1 keV dark matter particle, 
and we take y 2 = 12.939 as for thermal fermions (or 
sterile neutrinos produced via the Dodelson- Widrow 13011 
nonresonant mixing mechanism) and ga — 10.75 [lH, l40j| . 
This gives fc fs (r oq ) ~ 10/iMpc -1 . The comoving free- 
streaming distance, traveled by particles since t eq because 
of thermal velocities, 



Af s (t) 



_ch 

a(t>) 



(18) 



converges at late times to [3£ 



V6 



eq 



kfs (jeq ) 



0.24/i" 1 Mpc (19) 



This gives the scale below which the velocity dispersion 
of the warm dark matter has a strong effect. 



III. RESULTS 
A. Matter density power spectrum 

Since the equations of motion (fTTj) - (fT2"j) have the same 
form as those studied in [55| in the context of modified 
gravity models (but with a different kernel e(fc, r)) we use 
the same methods as in (55| for our numerical computa- 
tions. We refer the reader to [55| for a description of our 
analytical methods. 

Because of the explicit dependence on k introduced in 
the left-hand side of the Euler equation (fT2]l by the fac- 
tor e(fc, t), the linear growing and decaying modes of the 
density contrast now depend on k. As explained above, 
they also satisfy Eq.(|4]), with 5 = 0, where the depen- 

2 J 



dence on k is explicit. This pressurelike term, —k c 
slows down the growth of density perturbations at high 
k. For k > fcf s , which plays the role of a Jeans wave 
number, density perturbations would no longer grow but 
oscillate. We show in the lower panel of Fig. [1] the lin- 
ear growing mode D + (k,r) as a function of the wave 
number. (We normalize all linear growing modes to the 
CDM mode at the initial redshift Zi.) We clearly see 
the decrease of the growing mode above k ~ 1/Af s ~ Ah 
Mpc -1 . For comparison, we also plot in the upper panel 
the ratio P L , W dm / Pl ,cdm = T(k) 2 of the linear WDM 
to CDM power spectra, from Eq. (fT"5l) . All curves deviate 
from the CDM prediction at about the same wave num- 
ber, but, as expected, the damping of the linear power 
spectrum is stronger than the damping of the late-time 
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FIG. 1: Upper panel: ratio T{k) 2 — Pl,wdm/Pl,cdm of the 
linear WDM and CDM power spectra at z — 0, from Ea. (Tl5l) . 
Lower panel: linear growing mode D+(k, r) normalized to the 
scale factor a(r) as a function of the wave number. We show 
our results at z = and z = 3 for the CDM case (black dashed 
lines) and the WDM case (red solid lines) with m = 1 keV. 



linear growing mode. Indeed, the damping of the lin- 
ear power spectrum shown in the upper panel is mostly 
due to early-time effects, when the dark matter particles 
were still relativistic. In contrast, by definition of our 
initial conditions at = 100, the damping found in the 
lower panel is a late-time effect at z < Zi, due to the 
small nonzero velocity dispersion. This effect declines 
with time as the comoving wave number fcf s (r) grows as 
\fa from Eq.([6|). However, it is not zero and we will es- 
timate in the following the magnitude of this late-time 
effect. 

Next, we consider the nonlinear density power spec- 
trum in the perturbative regime in Figs. [5] and [3] We 
compare the results associated with the reference ACDM 
scenario, the WDM scenario with c s — (i.e., the only 
difference from CDM arises from the initial power (fTS]) ), 
and the WDM scenario with c s ^ (i.e., the differ- 
ence from CDM arises from both the initial power (fT5j) 
and the late-time velocity-dispersion) . To emphasize the 
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FIG. 2: Ratio of the power spectrum P(k) to a smooth ACDM linear power spectrum PLs(k) without baryonic oscillations, 
from [57| . We plot the linear power (lower group of dashed lines), the nonlinear one-loop "steepest descent" resummation 
(middle group of solid lines), and the "standard" 1-loop result (upper group of dotted lines). In each case, we show our results 
for the reference CDM scenario (slightly upper black lines), the WDM scenario with c s — (middle green lines) and with c s 7^ 
(slightly lower red lines). 
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FIG. 3: Relative deviation of the power spectrum from the CDM reference, in terms of the linear power spectra (APl/Pl, 
upper group of dashed lines), and of the nonlinear power spectra (AP/P, lower group of solid lines) obtained from the one-loop 
"steepest descent" resummation. In each case, we show our results for the WDM scenarios with c a = (slightly upper green 
lines) and with c s 7^ (slightly lower red lines). 



difference between various curves, we plot the ratios of 
our results by a common reference linear power spec- 
trum without baryonic oscillations, from [57] , in Fig. [5] 
We plot both the "standard" one-loop perturbative re- 
sult HU and the "steepest-descent resummation" (59r - 
l6lj , which agrees with the standard result up to one loop 
and contains a partial resummation of higher-order terms 

We recover the suppression of the nonlinear power 
spectrum due to the high-fc cutoff (TIB)) 46]. As in the 
usual CDM case, we can see that the nonlinear dynamics 
amplifies the power spectrum while erasing some of the 
baryonic oscillations. The difference between the "stan- 
dard" and the "resummed" perturbative predictions is 
similar to the one obtained in the CDM scenario, and it 
is larger than the difference between the CDM and WDM 



results. This means that on these scales, the standard 
perturbation theory is not accurate enough to describe 
the deviations between the CDM and WDM predictions 
(for to > IkeV), which are on the order of 1% as seen in 
Fig. [3] 

On the other hand, we can see that nonlinear contri- 
butions amplify these deviations, as compared with the 
linear power spectra. This is more clearly seen in Fig. [31 
where we plot the relative deviations from the reference 
linear and nonlinear CDM power spectra. The compari- 
son of the curves obtained with c s — and c s 7^ shows 
that most of the damping with respect to the CDM case 
is due to the cutoff (fi"6l) of the initial power spectrum. 
The late-time velocity dispersion only slightly amplifies 
the damping on these scales. Thus, for practical pur- 
poses, this late-time effect may be neglected. This jus- 
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FIG. 4: Halo mass functions for the reference CDM scenario (black dotted lines), the WDM scenarios with c s = (green dashed 
lines), and with c s 7^ (red solid lines). 



tifi.es the use of N-body simulations, initially built for 
CDM, to study gravitational clustering on these scales 
[ii], Hg]. However, for m > IkeV, the difference from the 
CDM power is not larger than 1.5% on these scales. This 
means that, as expected, most of the constraints that can 
be set from observations on the WDM scenario arise from 
smaller scales. However, because these smaller scales are 
also more difficult to predict (since they are deeper in 
the nonlinear or nonperturbative regime), it is interest- 
ing to check the signal that can be expected on the larger 
perturbative scales shown in Figs. [5] and which are bet- 
ter controlled since they can be described by systematic 
perturbative schemes. 

B. Halo mass function 

Another key statistic of large-scale structures is the 
mass function of collapsed halos. As in (5|| we define 
halos by a nonlinear density contrast of 200. Then, we 
obtain the halo mass function from the spherical dynam- 
ics, using the Press-Schechter scaling variable v [62l ]. 

with 

T~ x (200) 
v= q y \ (21) 

and the scaling function f{y) from [63| 

f(v) = 0.502 [(O.61/) 2 5 + (0.62z/) 05 ] e~ u2/2 , (22) 

which has been fitted to ACDM numerical simulations. 
This ensures that the halo mass function is always nor- 
malized to unity and obeys the large-mass tail n{M) ~ 
e~ u I 2 for any spherical-collapse mapping T q . We obtain 
the spherical dynamics associated with the equations of 



motion (fTTj) - (fT2")) as in (5f|. This provides the spherical- 
collapse mapping, S^q S x = J- q {5L q ), from the linear 
density contrast on the Lagrangian radius q to the non- 
linear density contrast on the Eulerian radius x. 

We show our results in Fig.|4j The low-mass tail should 
be considered with caution because its exponent may de- 
pend on the shape of the linear power spectrum and be 
different from the CDM case. Thus, numerical simula- 
tions suggest that a simple recipe of the form (l20l) , which 
involves a scaling function f{v) fitted to CDM simula- 
tions, overestimates the low-mass tail in WDM scenarios 
[44j . In contrast, the large-mass tail is better controlled 
because it is governed by spherically symmetric saddle 
points (and it does not involve the multiple mergers that 
affect the low-mass tail). 

As is well known, the WDM scenario leads to a much 
smaller halo mass function at low masses than in the 
CDM case, because of the lack of power on small scales 
[Til [44l [46| . For our purposes, Fig. 3] shows that the 
impact of the late-time velocity dispersion on the halo 
mass function at z < 5 can be neglected for m > IkeV. 
Indeed, the difference between the two WDM curves as- 
sociated with either c s ^ or c s = at z < Zi is much 
smaller than the deviation from the CDM reference. It 
is also smaller than the accuracy of halo mass functions 
that can be obtained from phcnomcnological models or 
numerical simulations. Again, this means that standard 
N-body codes can be used to predict the halo mass func- 
tion for WDM scenarios (with m > IkeV). 

This agrees with results from [g4| , which also model the 
velocity dispersion as an effective pressure term and com- 
pute the delay of halo collapse by this hydrodynamical- 
like pressure for a spherical dynamics. They also find 
that the initial cutoff of the linear power spectrum plays 
a greater role than halo dynamics, but the latter becomes 
significant at very low mass and high redshift (the late- 
time velocity dispersion has a very small effect below 
z < 40). 

Figure 2] can be compared with Fig. 7 of [65[ , which 
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FIG. 5: Probability distribution V(5 X ) of the matter density contrast within spheres of radius x, for three radii. We show our 
results at z = 3 for the reference CDM scenario (black dotted lines), the WDM scenarios with c s = (green dashed lines), and 
with c s 7^ (red solid lines). 



estimated the WDM halo mass function in the excursion- 
set formalism. These authors also find that the effect of 
the velocity dispersion, if it is estimated from its impact 
on the critical overdensity, is very small on the relevant 
mass scales (so that the cutoff of the halo mass func- 
tion at low mass is set by the initial cutoff of the power 
spectrum (|16|)). However, they find that if one modifies 
in addition the scaling function f{v), to take into ac- 
count the dependence of the mass function on the global 
shape of the function 5 C (M) (and not only on its value 
at the mass M of interest) , the low-mass cutoff becomes 
sharper. This second effect is not included in our paper 
(it requires additional modeling, such as the excursion set 
formalism, which needs to be checked against numerical 
simulations in this regime). 



C. Probability distribution of the density contrast 

Finally, we consider the probability distribution V(5 X ) 
of the density contrast within spheres of radius x. As in 
[5o| , we use a steepest-descent approach to obtain V{8 X ) 
from the spherical dynamics (gH, [66| . This is valid in the 
mildly nonlinear regime. We plot our results in Fig. [5] 
at z — 3, on scales that correspond to Lyman-a clouds. 
Indeed, depending on the details of the models Lyman- 
a clouds are associated with scales from ~ 10/i _1 kpc to 
~ l/i _1 Mpc (from the small Lyman-a forest clouds to 
damped systems) [67l - l69| . 

We recover the characteristic asymmetry induced by 
the nonlinear gravitational dynamics, with a shift of the 
peak toward underdensities (most of the volume is under- 
dense) , a very sharp low-density cutoff (6 X > — 1 since the 
matter density is always positive), and an extended high- 
density tail (most of the mass is within overdensities). On 
large scales, x > l/i _1 Mpc, all curves are very close and 
the deviations between the CDM and WDM scenarios are 
small. On smaller scales, the lack of small-scale power in 



the WDM scenario leads to a less advanced stage of the 
nonlinear evolution: the peak shifts closer to the mean 
(S x ) = 0, and the tails are sharper. The late-time velocity 
dispersion even further impedes the nonlinear evolution 
and makes the large density tail sharper. This suggests 
that accurate measures of the probability distribution of 
the flux decrement of distant quasars, due to Lyman-a 
absorption lines, which is closely related to the proba- 
bility distribution of the matter density on these scales 
[oTI l70j , could be sensitive to this late-time velocity dis- 
persion. Thus, numerical simulations that do not include 
this effect are likely to underestimate somewhat the dif- 
ference between the CDM and WDM scenarios with re- 
spect to Lyman-a absorption lines. 

It is not surprising that these statistics are more sen- 
sitive to the late-time velocity dispersion than the quan- 
tities studied in previous figures (the power spectrum on 
large perturbative scales and the halo mass function). 
Indeed, it is well-known that, because they probe rela- 
tively small scales, Lyman-a clouds are a sensitive probe 
of WDM scenarios. Most works focus on the decrease of 
the flux power spectrum due to the high-fc damping of the 
WDM power spectrum set by the relativistic free stream- 
ing jl3 - [2fj| . Figure [5] shows that the late-time velocity 
dispersion also has a non-negligible effect. 



D. Range of validity of our approximations 

Here we briefly comment on the validity of our method 
for the results described in the previous sections. Two 
effects are not fully included in our treatment of the dy- 
namics: 

a) nonlinearities beyond one-loop order for the compu- 
tation of the power spectrum, 

b) effects associated with the velocity dispersion that 
go beyond an effective pressure. 

The errors associated with a) can be estimated in Fig. [5] 
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FIG. 6: Left panel: relative deviation of the power spectrum from the CDM reference, as in Fig. using an initial redshift 
Zi = 50 (solid lines) or z% = 100 (dashed lines). Middle panel: ratio of the WDM halo mass function to the CDM halo mass 
function, with Zi = 50 (solid lines) or z, = 100 (dashed lines). Right panel: probability distribution of the matter density 
contrast, as in Fig. [5] for the CDM case, and for the WDM case with Zi — 50 (solid lines) and z% = 100 (dashed lines). In each 
panel, the green curves (which correspond to slightly higher P(k) and higher V(S X ) at S x > 2) are obtained for c s = and the 
lower red curves are for c s / 0. 



by the comparison between the standard perturbation 
theory and the resummation scheme. Even though the 
latter is systematically more accurate, the difference be- 
tween the two results gives an estimate of the error of the 
theoretical predictions on these scales. Then, as already 
noticed, one can safely conclude that, in this regime, the 
effect of the late-time velocity dispersion is on the order 
or smaller than the accuracy of standard perturbation 
theory. 

The errors associated with b) can be estimated from 
the difference between the cases c s — and c s > 0. As 
explained in Sec. Ill Al at late times, the main effect of 
the velocity dispersion should be well described by the 
"pressure" term in Eqs.(TJ]) and ©. Then, as long as 
the deviation between the two cases c s = and c s > 
is small, one expects that we obtain the correct order 
of magnitude. In the regime where the deviation would 
be large (at very high redshift), higher-order corrections 
are expected to come into play. Therefore, our method 
should be sufficient for our purposes, in the regimes stud- 
ied here(77j. 



IV. IMPACT OF LOWER INITIAL REDSHIFT 

In the previous figures, we set the initial conditions 
at redshift z, = 100. In this section, we investigate the 
impact of using a lower initial redshift, Zj = 50. This 
should make the two WDM results, with c s = and 
c s 0, closer to each other, since the velocity disper- 
sion decreases with time. However, because it also sets 
the nonlinear contributions to zero at Zi = 50 instead of 
Zi = 100, this also further underestimates gravitational 
clustering. 

In the standard ACDM scenario, we usually take the 



limit Zi — > oo within analytical approaches because one 
is only interested in accurate predictions [Hll 0, 
Moreover, numerical simulations often use second-order 
Lagrangian perturbation theory for their initial condi- 
tions to decrease the sensitivity to the initial redshift 
[73l ]. However, WDM simulations often use linear the- 
ory to set up their initial conditions. Indeed, a second- 
order implementation should, in principle, make use of 
a second-order analysis of the Vlasov equation at early 
times, while most works use the results obtained from 
the linearized Vlasov equation. On the other hand, the 
choice of the initial redshift is not so obvious in WDM 
scenarios, because a high initial redshift Zi can lead to 
spurious effects due to inadequate modeling of the large 
velocity dispersion [52j], whereas a low initial redshift al- 
leviates this problem but can lead to an underestimate 
of gravitational clustering. 

The dependence on the initial redshift Zi is easily in- 
cluded within our analytic approach as follows. In the 
perturbative framework, used for the large-scale power 
spectrum shown in Figs.[l]r21 the integrals over time asso- 
ciated with the one-loop and higher-order contributions 
run from Zi down to the redshift z of interest. In the 
spherical dynamics, used for the mass function and the 
density probability distribution shown in Figs. 0] and [SJ 
the equations of motion are also integrated from the ini- 
tial redshift z*. This allows us to obtain the dependence 
on zi of the large-scale structures built at a given redshift 
z, as in numerical simulations initialized at linear order 
at this redshift z%. We show our results for the density 
power spectrum, the halo mass function, and the density 
probability distribution function, comparing the choices 
Zi = 50 and z, = 100, in Fig-EQ. 

The left panel shows that on perturbative scales, k < 
0.9/iMpc" 1 at z = 5, using an initial redshift of z^ = 50 
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leads to an underestimation of the power spectrum on 
the order of 1%. Moreover, this is larger than the de- 
crease of the power spectrum due to the change from the 
CDM to the WDM power spectrum. On the other hand, 
the difference between the c s = and c s ^ results 
is very small. The same behavior is found for the halo 
mass function, as shown by the middle panel, with an 
underestimation of the large-mass tail because of the low 
initialization redshift z\ that is on the order of or larger 
than the true decrease due to the WDM scenario. The 
low-mass tail is less affected by the value of the initial 
redshift since it corresponds to moderate density fluctu- 
ations. In contrast, the large-mass tail corresponds to 
rare events, which amplifies the sensitivity to the initial 
conditions. This also explains why we find in the right 
panel that using a low initial redshift Zi = 50 does not sig- 
nificantly change the probability distribution of the den- 
sity contrast at z = 3, for moderate density fluctuations. 
Moreover, the difference between the c s — and c s ^ 
results remains similar to the one obtained with Zi = 100. 
Therefore, starting at z% = 50 instead of Z{ = 100 de- 
grades the accuracy of measures of gravitational cluster- 
ing at low redshifts, since it significantly underestimates 
the nonlinearities (as measured by the large-scale power 
spectrum or the halo mass function) and contaminates 
the signal associated with the WDM high-fc cutoff (on 
large perturbative scales or on the large-mass tail of the 
halo mass function). Moreover, using Zi = 50 does not 
help to reduce the effect of the late-time velocity disper- 
sion on the shape of V(6 X ) on Lyman-a scales. Thus, 
numerical simulations should use a high initial redshift, 
Zi > 100, rather than a low value, 2, < 50. 

V. CONCLUSION 

Using an effective Euler equation, that agrees with the 
Vlasov equation at the linear level (except for subdom- 
inant memory terms), we have estimated the impact of 
a late-time WDM velocity dispersion on the formation 
of large-scale structures. We have only considered the 
"cosmic web" , that is, large perturbative scales, moder- 
ate density fluctuations, and the number counts of virial- 
ized halos, which can be studied with analytic tools. We 
have focused on the case of a IkeV dark matter particle, 
which is representative of current WDM scenarios (lower 
masses are excluded by observations, such as Lyman-a 
forest data, while higher masses become indistinguish- 
able from the CDM limit). 

We find that on perturbative scales, the deviation of 
the density power spectrum from the CDM case is only 



on the order of 1%, at z < 5, even though it is slightly 
amplified by the nonlinear dynamics. This is below the 
accuracy of the standard perturbative expansion and re- 
quires efficient perturbative schemes. On the other hand, 
the effects of the late-time velocity dispersion are negli- 
gible over most of the perturbative range at z < 5 (so 
that one could use the same perturbative approaches de- 
vised for the CDM case). We also find that the late- 
time velocity dispersion has a negligible impact on the 
halo mass function at z < 5 (in any case, below the 
10% accuracy that can be guaranteed by simulations), 
although at very low mass and high redshift, the cutoff 
may become sharper [64l.l65l|. On the other hand, it has a 
non-negligible effect on the probability distribution of the 
density contrast on scales x < 0.1/i _1 Mpc at z = 3. This 
means it should have some impact on the probability dis- 
tribution of the Lyman-a flux decrement, measured on 
the spectra of distant quasars. 

Finally, we note that numerical simulations should use 
a high initial redshift, Zi > 100, rather than a low value, 
Zi < 50. Indeed, such a low initial redshift can lead to 
a significant underestimation of the power spectrum on 
perturbative scales and of the large-mass tail of the halo 
mass function, which is larger than the true signal asso- 
ciated with the WDM scenario (but of course, on smaller 
scales and on the low mass tail of the mass function, one 
is again dominated by the actual WDM signal). A low 
initial redshift does not help much either to reduce the ef- 
fect of the late-time velocity dispersion for the probability 
distribution of the density contrast on scales associated 
with Lyman-a clouds. 

To go beyond the effective hydrodynamical equations 
used in this work, one should use the nonlinear Vlasov 
equation itself. However, this is a difficult task because 
of the additional velocity coordinates, which makes nu- 
merical implementations significantly heavier already for 
the CDM scenario 0, [75[ . An alternative would be to 
extend the fluid approximation to higher orders (50j . by 
including equations of motion for the velocity moments 
of the Vlasov equation up to some higher order n > 3. 
However, because most of the WDM signal arises from 
nonperturbative scales, such a task may not be very re- 
warding, unless one builds methods that can be applied 
to the Lyman-a forest clouds, for instance. 
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